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Abstract 

A review on the two-pion Hanbury-Brown-Twiss (HBT) interferometry in 
the granular source model of quark-gluon plasma droplets is presented. The 
characteristic quantities of the granular source extracted by imaging analy- 
sis are presented and compared with the HBT radii obtained by the usual 
Gaussian formula fit. The signals of granular sources are presented. 



1 Introduction 

The main purpose of relativistic heavy ion collisions is the study of the quark- 
gluon plasma (QGP) formed in the early stage of the collisions. As is well 
known, the system produced in the collisions may thermalize and reach local 
equilibrium in a very short time To- The following evolution can be described 
by hydrodynamics. Hydrodynamics provides a direct link between the early 
state and final observables. It has been extensively used in relativistic heavy 
ion collisions. However, the usual hydrodynamic calculations cannot explain 
the RHIC HBT data, R out /R side « 1, the so-called RHIC HBT puzzle pQ. 

In reference [2] a granular source model of QGP droplets was put forth for 
the HBT puzzle. In references [31 H] the granular source model was improved 
to explain the RHIC HBT data pQ. In reference [5J, the fluctuations of the 
single-event two-pion correlation functions of granular sources were discussed. 
In references [6j [7], some observables for granular sources were put forth 
and analyzed for the source inhomogeneity in the heavy ion collisions at 
RHIC and Large Hadron Collider (LHC) energies. The imaging characteristic 
quantities and signals of the granular source of QGP droplets were presented 
in references [U [5] . This paper will give a review on the important ingredients 
and recent progresses of the granular source interferometry. 

2 Granular source model of QGP droplets 

The granular structure of QGP droplets may occur in the first-order phase 
transition from the QGP to hadronic gas, which was proposed by E. Witten 
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in 1984 [9]. In references [TUl EEE], the two- and multi-pion HBT correlation 
functions of granular sources were investigated. 

In two-pion interferometry, there is the relation among the HBT radii 
R out and R S id e , the transverse velocity of pion pair vt, and source lifetime r, 
Rl ut ~ -Rfide + v r 7 " 2 [12] • I n reference [2], it was noticed that the lifetime r for 
a uniform hydrodynamic evolving source scales with the source initial size. 
A smaller (or larger) initial source radius has a smaller (or larger) -R S id e , and 
smaller (or larger) source lifetime r. So, there is always R out > -Rside for the 
uniform hydrodynamic sources. However, for a granular source with small 
QGP droplets distributed in a large region, the source HBT radius R out is 
approximately equal to R S id e - It is because that for the granular source the 
HBT radius -R S id e (proportional to the QGP droplets distribution) is much 
larger than the source lifetime r (proportional to the droplet radius). 

In relativistic heavy ion collisions, the system initial transverse energy 
density is highly fluctuating on event-by-event basis. The large initial fluctu- 
ations together with the effects of violent expansion and surface tension may 
lead to formation of granular droplets [3J [32] • On the other hand, the bulk 
viscosity of the QGP may increase rapidly near the QCD transition, which 
also leads to formation of the QGP droplets [14]. 

It is assumed in the improved granular source model [4] that the system 
produced in central relativistic heavy ion collisions fragments and forms the 
QGP droplets at a time t after r . The droplets distribute initially within a 
short cylinder along the beam direction (z direction) with the probabilities 
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1 - exp (-pl/AKl) 0(K±-po), (1) 



dP 

-T 1 = 0{y m - 1 2/o | ) , zq = t tanh?/o, (2) 

where po an d z are the initial transverse and longitudinal coordinates of the 
droplet center, y is the initial rapidity of the droplet, TZ±, A7Z±, and y m are 
the radius of the cylinder, the shell parameter of the initial distribution, and 
the central rapidity limitation. On the basis of the Bjorken picture [15], the 
velocity of the droplet is assumed as 

Vd± = ^(^/V 1 ~ v dz> V dz = z o/t , (3) 

where a? and bx are the magnitude and exponential power parameters. 

The evolution of the system after the fragmentation is the superposition 
of all the evolutions of the individual droplets, each of them is described by 
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relativistic hydrodynamics with a cross-over equation of state of the entropy 
density [16J. The initial radius r' Q of the droplets in local frame is assumed 
having a Gaussian distribution with standard deviation a. To include the 
pions emitted directly at hadronization and decayed from resonances later, 
the pion freeze-out temperature is taken in a wide region with the probability 

dP f 

cx / dir exp[-(T - T h )j AT dir ] + (1 - / dir ) exp[-(T - T h )/AT dec ] , (4) 

where /dir is a fraction parameter for the direct emission, is the tem- 
perature of complete hadronization, AT dir and AT dec describe the widths of 
temperature for the direct and decayed pion emissions. 

The selections for the model parameter values can refer to [3]. 



3 Granular source imaging 

The imaging technique introduced by Brown, Danielewicz, and Pratt [T7j 
is a model-independent way to extract the two-pion source function S(r), 
the probability for emitting a pion pair with spatial separation r in the pair 
center-of-mass system (PCMS), from the HBT correlation function. 
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Fig. 1. Three-dimension source 
functions of the granular source. 
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Fig. 2. (a)-(d) Imaging results of granu- 
lar source, (e)-(f) HBT radii of granular 
source and RHIC experiments pQ. 



Figure 1 (a), (b), and (c) show the source functions of the granular source 
in x (out), y (side), and z (long) directions. Here k T is the transverse mo- 
mentum of pion pair in the longitudinally comoving system (LCMS). One 
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can see that the source functions in out and long directions have long tails. 
The reason is the source expansion which boosts the pion pair in out and 
long directions. The source function width in long direction is smaller for 
larger kx because the average longitudinal momentum of the pairs is smaller 
for larger kr- 

Once obtain the source functions one can introduce the quantities of r 
moments to describe the source geometry numerically. They are model- 
independent quantities. Figure 2 (a)-(c) exhibit the quantities, Ri =y/n(ri) /2, 
(i = x,y,z), of the first-order moment (r,) for the granular source. Ri de- 
scribes the source size in i direction and normalized to the Gaussian radius 
for one-dimension Gaussian source [I]. For comparing, the usual Gaussian 
fitted HBT radii of the granular source and RHIC experiments are shown in 
Fig. 2 (d)-(f). In Fig. 2 (a), the symbols v denote the results of •y^ 1 R x , 
where 7y is the Lorentz contracted factor of LCMS to PCMS. One can 
see that the imaging results of Ri are consistent with the HBT radii after 
considering the Lorentz contraction. 



4 Granular source signals 

For granular sources, the single-event HBT correlation functions exhibit large 
event-by-event fluctuations [SI El El- In order to observe the event-by-event 
fluctuations, we introduce the quantity / of the relative fluctuation of the 
single-event correlation function to mixed-event correlation function, with its 
error-inverses weight [6] 
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Fig. 3. The distributions dN/df for 
40 events with FIC and SIC. 
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Fig. 4. The root-mean-square of 
/ as a function of N n7T . 
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Figure 4 shows the distributions of the / for the variables of transverse rela- 
tive momentum gtrans and relative momentum q of the pion pairs for the 40 
events generated by the smoothed particle hydrodynamics US], with the 
impact parameter 6 = 5 fm. It can be seen that for the pion pair number 
N nn = 5 x 10 6 , the distributions for the fluctuating initial conditions (FIC) 
are much wider than those for the smoothed initial conditions (SIC) both for 
(ftrans and q. Even for N nn = 5 x 10 5 , the widths for FIC are visibly larger 
than those for SIC. Figure 5 shows the root-mean-square (RMS) of /, as a 
function of N n7T for the events. It can been seen that the values of / rms rapidly 
increase with N nn for FIC because the errors in Eq. §5§ decrease with N nn . 
For SIC the values of / rms are almost independent from N nn , because both 
the differences and their errors in Eq. <^ decrease with N n7T in the case. At 
LHC energy the event multiplicity of identical pions is about two thousands 
and the order of N nn will be 10 6 . In this case, the distributions of / and its 
RMS will provide observable signals for the source inhomogeneity. 



5 Discussion and conclusion 

In relativistic heavy ion collisions, the fluctuating initial density distribution, 
violent expansion, surface tension, and viscosity of the QGP may lead to 
formation of the granular sources. For the granular sources, the single-event 
HBT correlation functions have large event-by-event fluctuations. Because 
of data statistics the experimental HBT analyses are performed with the 
mixed-event HBT correlation functions, and the fluctuations are smoothed 
out after the event average. However, the short lifetime of the granular source 
may survive and lead to the HBT puzzle. The model-independent quantities 
Ri obtained by imaging analysis for the granular sources are consistent with 
the usual HBT radii of the granular sources and RHIC experiments. At the 
coming LHC heavy ion collisions, it is hopefully to observe the signals of the 
granular sources. 
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